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Abstract. The global missing baryon problem — that the sum of observed baryons falls short 
of the number expected form BBN — is well known. In addition to this, there is also a local 
missing baryon problem that applies to individual dark matter halos. This halo by halo missing 
baryon problem is such that the observed mass fraction of baryons in individual galaxies falls 
t ' short of the cosmic baryon fraction. This deficit is a strong function of circular velocity. I give 
' an empirical estimate of this function, and note the presence of a critical scale of ~ 900 kms~^ 
; therein. I also briefly review Q.f, from BBN, highlighting the persistent tension between lithium 
and the CMB, and discuss some pros and cons of individual galaxies and clusters of galaxies as 
potential reservoirs of dark baryons. 
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Q-i' 1. Introduction 
I 

O Modern cosmology suffers a hierarchy of missing mass problems. Most famously, there 

appears to be more gravitating mass than can be accounted for with baryons while 
maintaining consistency with Big Bang Nucleosynthesis (BBN): 



Vlb < n„,. (1.1) 



J> ■ This, together with the need to grow large scale structure, leads to the inference of 
lO ' non-baryonic cold dark matter (CDM). 

In addition to the dynamical missing mass problem, there is also the missing baryon 
problem (e.g.. lPersic fc Salucci 1992l|Fukugita et al. 1998p . Our inventory of known baryons 
in the local universe, summing over all observed stars, gas, etc., comes up short of the 
total anticipated from BBN: 



o 

. n^ + ng + ... <nb. (1.2) 

O 



For example. Bell et al. (2003) estimate that the sum of stars and cold gas is only ^ 3% 
of fif,. While there now seems to be a good chance that many of the missing baryons 
^> \ are in the form of highly ionized gas in the warm-hot intergalactic medium (the WHIM 
^ ' — see Mathur, these proceedings), we are still far from being able to give a confident 
5^ \ accounting of where all the baryons reside. Indeed, there could be multiple distinct reser- 
voirs in addition to the WHIM, each comparable to the mass in stars, within the current 
uncertainties. 

Here I highlight the halo by halo missing baryon problem. In addition to the global 

missing baryon problem, there seems to be a mismatch between the mass of baryons in 
individual galaxies and the mass of their host halos. For the amount of dark matter we 
infer dynamically, the mass in detected baryons falls short of the cosmic baryon fraction: 

TUd < fb- (1.3) 
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Figure 1. The baryon density from various measurements over the past decade, as tabulated by 
McGaugh (2004), with a few recent updates (|Ch arbonnel fc Primas 2005' 'Asplund et al. 2006"1 
[Bonifacio It al. 20071 lO'Meara et al. 20061 [Peimbcrt et al. 2007, Spcrgcl et al. 2007J). Con- 
straints from each independent method are noted by different symbols. BBN was well established 
long before the start date of this graph; previous work is represented by early compilations. This 
plot assumes Ho = 72 kms~^ Mpc^^; the first point is the famous Qbh^ ~ 0.0125 value of Walker 
et al. (1991). Note that no measurement of any isotope ever suggested a value Q,bh^ > 0.02 (hor- 
izontal dotted line) prior to the appearance of relevant CMB data (vertical dotted line). With 
the exception of the case of no CDM (McGaugh 1999, 2004), no CMB fit has ever suggested 
Qbh^ < 0.02. Measured values of some isotopes seem to have drifted upwards towards the CMB 
values since 2000, a worrisome trend given the historical tendency for cosmological measurements 
to track one another. 



Here 
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(1.4) 



is the mass fraction of detected baryons in a given object, and the cosmic baryon frac- 
tion is /b = 0.17 dSpergel et al. 2007D . Unfortunately, the virial mass of any given galaxy 
is a rather notional quantity, and there are many possible mappings between the ob- 
served baryonic content and the inferred dark mass. Nevertheless, there is a pronounced 
regularity to the data that does not tumble naturally out of our ACDM theory. 



2. Big Bang Nucleosynthesis 

Before proceeding with a discussion of the missing baryon problem in individual dark 
matter halos, it is worth reviewing SI;,. It is commonly stated that we have entered the 
era of precision cosmology. If this is true, then there is a tension between fif, — 0.03 and 
0.04 that should concern us. 

One of the great successes of cosmology is BBN. The abundance of the light elements 
depends on a single parameter, the baryon-to-photon ratio, that maps directly to the 
baryon density. Independent measurements of ■^H, "^He, and ^Li all indicate nearly the 
same value for fib- In addition, modern observations of the cosmic microwave background 
(CMB) provide another similar constraint on the baryon density. 
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Table 1. fib by method*. 



W9lt 


2H 


^He 




CMB 


W3t 


M04* 


0.024 
a 0.005 


0.038 
0.005 


0.021 
0.008 


0.028 
0.005 


0.046 
0.007 


0.0430 
0.0014 


0.033 
0.006 



'Assumes Hq = 72 kms ^ Mpc ^ . 
Compilation of Walker et al. (1991). 
*WMAP-only ACDM fit from A website. 
*CMB without CDM — see McGaugh (2004). 



The consistency of the independent BBN constraints is impressive, and I have little 
doubt that the basic picture is correct. However, the agreement is not perfect. There 
seems to be a modest dichotomy between each of the independent methods (Fig. [T]), 
with each preferring a slightly different value of Table [T] shows a comparison of 
baryon densities by method. Values are the median for each case from the compilation 
of |McGaugh 2004| assuming = 72 kms-iMpc"^ The uncertainty is taken from the 
variance of different published measurements. This tends to overstate the error claimed by 
individual determinations. For this uncertainty, the light elements are broadly consistent. 
If we look only at the latest WMAP result ( Spergel et al. 2007| ), its tiny uncertainty 
makes it difficult to reconcile with the helium and lithium. I briefly review each method 
below. 

Helium: Measured in very low metallicity extragalactic HII regions (e.g., Kuzio de Na- 
ray et al. 2004), helium gives the lowest value of Slf,. Though consistently in the vicinity 
oi Yp = 0.24, the weak dependence of Yp on 17f, makes it difficult to obtain a precise 
estimate from helium. Moreover, independent determinations vary by more than they 
should for the stated errors. This is likely due to the many subtle systematic effects 
that afflict helium abundance determinations at the ~ 1% level (jOlive fc Skillman 2004p . 
Thus, while the broad consistency of helium with other determinations is comforting, it 
seems unlikely that it can provide a strong constraint on the precise value of f2b. We 
should not take this as a free pass to simply ignore helium; all determinations of Yp up 
to the date of my review gave fib < 0.035, and I found it increasingly difflcult to reconcile 
any of the helium data with Q,b > 0.04 (but see lPeimbert et al. 2007p . 

Lithium: Measured in low metallicity stars (e.g., |Ryan et al. 2000 ), lithium gives the 
next lowest value for fib. Unlike helium, the ^Li abundance seems to be both well mea- 
sured and consistent between independent determinations. The Spite plateau is sharply 
defined by dozens of low metallicity stars. At fif, = 0.028, it is marginally consistent 
with deuterium but not consistent with the CMB. The general presumption seems there- 
fore to be that lithium needs to be 'fixed' (e.g.. lKorn et al. 20061 IPiau et al. 2006p . but 
there is nothing broken about it if the CMB value is overestimated or its uncertainty 
underestimated. 

Deuterium: Of the light elements, deuterium is the most sensitive to fib- Its mea- 
surement in Lyman absorption clouds in the spectra of high redshift quasars provides 
one of the most accurate constraints on BBN. Indeed, it came in at surprisingly low 
D/H, giving high Qbh'^ = 0.019 ± .0012 (|Tytler et al. 19991 ) relative to long-standing 
resuhs (e.g., ^th'^ = 0.0125 ±0.0025: IWalker et al. 1991|) . Unfortunatelv. there exist only 
6 accurate, published D/H values (fewer systems than there are points in Fig.[T]), one of 
which is discrepantly high. 

CMB: The microwave background provides another constraint on flf, that is indepen- 
dent of abundance measurements. It gives a value for Qb that is close to but persistently 
higher than all of the light elements. The power spectrum of the CMB is very sensitive 
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Figure 2. The effects of the baryon density and baryon fraction on the CMB. The points are 
the third year WMAP data (Spcrgcl et al. 2007). The various hnes are models with Q.b and fb 
set to fixed values (see text). The dotted (topmost) line is set to the lithium value {0.^ = 0.028) 
with fb — 0.12 (case 1). The dashed (middle) line is set to the deuterium value {Qb — 0.038) with 
fb ~ 0.17 (case 2). Case 3 is the solid (bottom) line with lithium fib and no CDM (McGaugh 
1999, 2004). The CMB is exquisitely sensitive to both the absolute baryon density and the 
baryon fraction, but also suffers degeneracies between these and other parameters (especially 
the tilt). 

to both the absolute density of baryons (ilf,) and the global baryon fraction {fb)- Fig. [5] 
illustrates this sensitivity with three cases: (1) a plausibly low baryon density, low baryon 
fraction universe; (2) a high baryon density, high baryon fraction universe, and (3) a low 
density, no-CDM universe. 

Case (1) is motivated by the lithium baryon density (Table [1]). For this value of fib, 
fb = 0.12 as observed in clusters ( [McCarthy et al. 2007[ see also Holder, these proceed- 
ings) if flm = 0.23, consistent with peculiar velocity measurements (fi™ = 0.22 ± 0.02: 
|Moyahee fc Tully 2005p . These three independent lines of evidence provide a self consis- 
tent and reasonable set of cosmological constraints. However, they badly overestimate 
the amplitude of the second acoustic peak. This problem is generic to all concordance 
ACDM models prior to the first measurement of the second peak in 2000. 

Case (2) utilizes the deuterium abundance, which gives the highest fib of all the BBN 
elements and the one most consistent with fits to the CMB. A baryon fraction fb — 0.17 
is adopted, giving a plausible fl„i = 0.26. Though rigged to be as consistent as possible 
with CMB fits without simply adopting them, this case still over-predicts the amplitude 
of the second peak, albeit by a lesser amount than case (1). 

Case (3) uses the lithium fib with zero CDM. Consideration of this model is motivated 
by MOND, but uses entirely conventional physics. It provides quite a good fit to the 
second peak. Indeed, this was the only successful a proiri prediction of its amplitude 
( [McGaugh 1999p . Here it merely illustrates the sensitivity to as well as fib- 

Like all models devoid of CDM, case (3) under-predicts the amplitude of the third 
peak, just as un-tilted ACDM models over-predict it (e.g., cases 1 and 2). This 'low' third 
peak leads to the inference of a tilt (n w 0.95), which in turn impacts the inferred baryon 
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density, somewhat reducing it relative to previous CMB determinations. In the context of 
MOND, the 'high' third peak impHes a driving term analogous to that provided by CDM. 
This occurs naturally in TeVeS-like relativistic MOND theories (|Skordis et al. 2006|) , but 
it is unclear how this might affect fli,. 

2.1. BBN summary 

BBN is undoubtedly one of the shining successes of modern cosmology. There is broad 
agreement between multiple independent lines of evidence to an accuracy rarely witnessed 
in cosmology. I have little doubt that the basic picture is correct. 

In detail, there are significant differences between independent methods. It is common 
to point to the great consistency between from BBN and the CMB as support for 
the ACDM picture. However, this is rather sloppy, being only really true for deuterium. 
Lithium and helium prefer noticeably lower values for the baryon density. 

The light elements themselves provide a self-consistent picture if we adopt the uncer- 
tainties listed in Table [TJ It is the tiny uncertainty in the CMB determination that is 
problematic. The CMB is indeed that sensitive (Fig.[2|), but does suffer from parameter 
degeneracies, in particular the tilt and the baryon fraction. The tension between the light 
elements and the CMB disappears if there is no CDM. 

It appears that the true baryon density is not yet as precisely known as might be hoped. 
Nevertheless, the missing baryon problem remains an issue at low redshift. Whether the 
true baryon density is 3% of critical or 4% depends on whether one gives more weight 
to light element measurements or CMB fits. Either way, fif, exceeds the well-determined 



inventory of local baryons (Fukugita et al. 1998) 



3. The baryon content of individual dark matter halos 

As well as trying to account for all the baryons in the universe, and the integrated 
amount in galaxies, clusters, and the IGM (see Frenk; Moore, these proceedings), it would 
be nice to have an accounting in individual systems. Each dark matter halo can, to a first 
approximation, be thought of as a microcosm of the whole. As such, one would naively 
expect each halo to have the same baryon fraction as the whole universe. On the scale 
of clusters of galaxies, this is approximately true. For individual galaxies, observations 
depart from this ideal in a way which we have yet to understand. 

3.1. The Milky Way 

Our own galaxy provides an illustration of the general problem. We have a rather com- 
plete inventory of local baryons that appears to be fairly consistent between independent 
estimates. Yet when we compare the baryonic mass to the total dynamical mass, it falls 
well short of the universal baryon fraction. 

For illustration, I adopt estimates of the baryonic mass from the recent work of 
Flynn et al. (2006)| Mf^ « 5 x 10^° M© and Mf^ « 10^° Mq. Most of the stel- 



lar mass is in the disk, though the spit between disk and bulge depends somewhat on 
details like the disk scale length. The gas mass is split between molecular and atomic 
components. Other baryonic components appear to be negligible. The total baryonic 
mass of the Milky Way is thus Mf ^ « 6 x 10^° Mq. 

The baryonic masses of satellite galaxies within the Milky Way halo add little to this 
sum. Their total is probably less than the uncertainty in the Milky Way mass, and 
certainly less than 10^" M©. It is thus hard to imagine that the baryonic mass associated 
with the Milky Way halo is any larger than 10^^ M©. 

The total mass of the Milky Way has been studied many times. The answer persistently 
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Figure 3. Rotation curves for two galaxies illustrating the circular velocity measured from 
the outer regions of extended 21 cm data (horizontal lines). The distribution of the baryonic 
component is also illustrated. For the purposes of this review, we are only concerned with the 
total mass in baryons and the characteristic circular velocity of the system. 



comes back in the vicinity of M. 



MW 



2 X 10^^ M0. There are many uncertainties in this 



estimate, the dominant one being whether the most distant tracer, Leo I, is considered 
to be bound to the Milky Way or not. ISakamoto et al. (2003)1 give 2.5 x 10^^ Mq with 



Leo I and 1.8 x 10^^ Mq without it. A practical lower limit is 10^^ Mq; half this is 
already encompassed within the radius of the LMC. 
Taking these numbers at face value. 



AlW 



6 X lO^" M, 
2 X 1012 M, 



0.03. 







(3.1) 



This falls well short of the universal baryon fraction, fi, — 0.17. Pushing the numbers 
to their extremes, we can place the limit m^^^ < 0.10. This is at least within hailing 
distance of fb — 0.12 from clusters, but still well short of the CMB value. 

This raises the question: where are all the baryons that should be associated with the 
Milky Way halo? Perhaps they have been expelled, though the Milky Way is generally 
thought to be too massive to have suffered much of this. Another possibility is that they 
lurk in the halo in some as yet undetected form, like a local WHIM. While the WHIM is a 
great candidate for many of the globally missing baryons (see Mathur, these proceedings), 
it appears to fall well short of what is needed to explain this deficit in the Milky Way 
dPildis &: McGaugh 1996| , providing perhaps ~ lO^" Mq (jCollins et al. 2005|) . 

The Milky Way is not unique in having < /;,. Similar discrepancies are observed in 
many systems ( [McGaugh fc de Blok 1998] ). Indeed, nid appears to deviate systematically 
from /f, by an amount which becomes increasingly severe in systems of decreasing circular 
velocity. 

3.2. The Baryonic Tully-Fisher Relation 

A basic tenet of CDM theory is a relation between halo mass and circular velocity at the 
virial radius (jSteinmetz &: Navarro 1999|) of the form 



(3.2) 



with the normalization depending on the cosmology (|Mo &: Mao 2004]) . We can measure 
neither of these quantities. What we can measure is the baryonic mass of a given galaxy 
and a characteristic circular velocity (Fig. [3]). 

Summing up the mass in stars and gas for a large sample of spiral galaxies, we obtain 
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Figure 4. The baryonic mass-circular velocity relation, spanning the range from the tiniest 
dwarfs (squares) through spiral galaxies (circles; McGaugh 2005) to clusters of galaxies (trian- 
gles; |Reiprich 2001 1 . The dashed line is the fit to the galaxy (circles only) data; the dotted line 
is what we expect in ACDM if all the baryons in each halo are identified. 



the Baryonic TuUy-Fishcr relation: 



Mb = 50K* 



(3.3) 



with Mb in solar masses and Vc in kms ^ (Fig. [4]). For an extensive discussion of the 
effects of various choices for the stellar mass-to-light ratio on the normalization, slope, and 



scatter in this relation, see McGaugh (2005), Equation 13.31 does well when extrapolated 
to predict the locations of systems with Vc < 50 kms^^, including both rotating gas rich 
dwarf Irr systems (the squares in Fig. [4]) and pressure supported, star dominated dwarf 
spheroidals (Mateo & Walker; Geha, private communications). 

Extrapolation of the Baryonic TuUy-Fisher relation to the cluster scale does less well, 
over-predicting the median cluster mass by a factor of ~ 2.5. Perhaps these are simply 
different systems, and such a vast extrapolation is unwarranted. On the other hand, 
coming this close with an extrapolation over several decades in mass is fairly remarkable. 

That clusters fall a little below the extrapolation of the Baryonic TuUy-Fisher relation 
may suggest that our census of baryons is not complete in these systems. If so, clusters 
might be a good place to look for some of the missing baryons. Clusters are not a large 
contributor to the total integrated baryonic mass of the universe, so even if more baryons 
are lurking there, they are not likely to solve the missing baryon problem. Placing them 
on the extrapolation of equation 13.31 would require {very roughly) 5% of fib- While not 
a large fraction of the total, this is comparable to all the stars in field galaxies. Most of 
the missing baryons could be in the WHIM and still leave room for a number of other 
distinct baryonic reservoirs like this. 

The CDM mass-velocity relation comes much closer to the cluster data than does the 
empirical Baryonic TuUy-Fisher relation. The difference between it and the median cluster 
value is the difference between fb = 0.17 and 0.12. This near agreement deteriorates for 
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Figure 5. The baryon content as a function of circular velocity. Points are galaxy an d cluster 
data, as per Fig. |4l The solid line is the reionization model from Grain et al. (200"7)| (see also 
[Mayer fc Moore 2004) ; the dashed hne is the preheating model of Mo & Mao (2004); the dash-dot 
line is a fit to the data (see text). 



smaller systems. The slope of the halo mass-virial velocity relation is shallower than the 
observed baryonic mass-circular velocity relation. Even bearing in mind the caveat that 
the halo virial and observed circular velocities may not be simply related, it is difficult 
to avoid the conclusion that the observed number of baryons is systematically less than 
the expected number as one goes down the mass function. 

If we take the halo mass-virial velocity relation as sacrosanct, then it is impossible to 
avoid the conclusion that varies with Vc (Fig. [5]). Empirically, this variation is well 
fit by 

rrid = fb tanh . (3.4) 

Here V| is a scale where nid — > fb', the dash-dot line shown in Fig. [5] uses equation 13.41 
with = 900 kms"^. This is the scale where the dashed and dotted lines cross in Fig. [5] 
which is where the data transition from looking like ACDM on larger scales to looking 
like MOND on smaller scales. Note that something odd is happening on intermediate 
scales, as groups in the Local Volume imply flm ~ 0.04 ([Karachentsev 2005,) . (Yes, that 
is ^Im, not fib-) 

In order to derive eauation l3.4l or plot Fig. [5l it is necessary to make some assumption 
in order to relate total halo mass to observed baryonic mass. For simplicity, I assume 
Vvir — Vc here. While I do not expect this to hold in detail, plausible adjustments (due, 
for example, to adiabatic compression) will only tweak the details and not alter the 
basic result. That equation 13.41 returns a value consistent with that for the Milky Way 
(eauation l3.1[) suggests that the simple assumption is not terrible. 

The variation of nid with Vc implies that there should be missing baryons associated 
with each and every galaxy halo (see also Zackrisson, these proceedings). Equation 13.41 
is merely a fit to the data. What should we expect? Recent work gives some idea: Mayer 
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& Moore (2004) and Grain et al. (2007) consider the effects of reionization while Mo & 
Mao (2004) consider the possibihty of preheating in analogy with clusters. Unfortunately, 
neither of these types of models matches the observed trend (Fig. O. In particular, there 
is no evidence in the data for the sudden transition expected around 20 kms"-'^ from 
reionization. 

In the context of ACDM, this merely means that there is more physics than reionization 
mediating between the primordial mass-velocity relation and the observed one. Feedback 
is frequently invoked in this context, but this word has come to mean many different 
things. It remains to be demonstrated whether some flavor of feedback can be a solution 
to all these problems. 

There are, in principle, many potential sources of scatter that should manifest them- 



selves in the Tully-Fisher relation (McGaugh & de Blok 1998), yet the observed rela- 



tion has remarkably little scatter. It is hard to imagine how nid can vary so much be- 
tween bright and faint galaxies with virtually zero scatter along the way. Somehow, a 
Vc ~ 50 kms~^ galaxy must 'know' it should have « 0.01 and never 0.04, while for 
a galaxy with Vc ~ 220 kms""'^ the situation is reversed. The small observed scatter only 
happens naturally if = /b ( |McGaugh et al. 2000 ). This implies that essentially all the 



baryons are accounted for on galaxy scales, making individual galaxies an unlikely place 
to look for the missing baryons. However, this contradicts a basic tenet of CDM, as the 
intrinsic mass-velocity relation differs from equation 13.21 if nid = fb- 



4. Conclusions 

I provide two sets of interpretations, one in the context of ACDM, the other from an 
empirical perspective. 

4.1. Conclusions in the context of ACDM 

• There is no signature in the data due to reionization in the vicinity of ^ 20 kms^^. 
Low mass galaxies smoothly follow the extrapolation of the Baryonic Tully-Fisher relation 
fit to galaxies with Vc > 50 kms~^. 

• The baryon content of halos scales with circular velocity as summarized by equa- 
tion O 

• Models invoking reionization or pre- heating fail to match the observed trend in m^. 

• Physics beyond reionization is required to explain the observations. 

In order to reconcile the data with AGDM, some physical effect must make the observed 
baryon fraction the observed function of circular velocity. Whatever mechanism is 
invoked must operate with remarkably small scatter, posing a fine-tuning problem. The 
physics of galaxy formation — the bridge connecting the primordial halos found in nu- 
merical simulations with real observed galaxies — remains as sound and secure as the 
Bridge of Death over the Gorge of Eternal Peril in Monty Python's Holy Grail. 

4.2. Empirical conclusions 

• The observed mass-velocity relation is steeper than predicted by CDM. 

• The small scatter in the Baryonic Tully-Fisher relation implies that substantially all 
the baryonic mass has been accounted for in individual galaxies. 

• The two preceding items contradict a basic tenet of CDM. 

• Extrapolation of the Baryonic Tully-Fisher relation to cluster scales suggests that 
the inventory of baryons in clusters may be incomplete. 

The last point implies that perhaps some of the missing baryons are hiding in clusters in 
some yet to be identified form. These need not be a large fraction of the missing baryons. 
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maybe 5%. There is precedence for such a situation. Nowadays it is widely accepted 
that the hot X-ray emitting intra-cluster gas is the dominant baryonic component of 
rich clusters. A couple of decades ago, stars were widely thought to be the dominant 
baryonic component. Though I am at a loss to suggest a plausible state in which baryons 
might hide (but see IBonamente et al. 20031) . this history reminds us that we should not 
be over-confident about the non-existence of matter we haven't yet detected. 
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